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1.  OBJECTIVES 

The  contractor’s  objective  in  the  first  year  is  to  develop  a  technique,  based  on 
minimum  entropy  deconvolution  (MED),  useful  for  discriminating  ripple-fired  explosions 
from  other  seismic  events. 

2.  PROGRESS 

An  interactive  X-Windows  program,  using  modules  developed  for  the  Nuclear  Moni¬ 
toring  Research  and  Decelopment  (NMRD)  initiative,  has  been  developed  to  analyze  the 
effects  of  MED  filters  on  explosion  data.  Using  this  program,  we  have  started  to  evaluate 
the  effects  of  the  use  of  various  parameters  in  the  problem:  filter  length,  damping 
coefficient,  data  window  position  and  length.  The  algorithm  for  generating  an  MED  filter 
is  based  on  an  iterative  search  for  the  local  minimum  of  a  specific  norm  of  the  windowed 
data.  The  program  permits  interactive  analysis  and  evaluation  of  the  evolution  of  the 
filter,  as  well  as  its  effect  on  the  data,  after  each  iteration  by  visual  monitoring  and  inter¬ 
cession.  We  have  found  that  the  "best"  filter  is  frequently  encountered  before  a  local 
minimum  is  reached. 

We  have  examined  signals  from  the  Soviet/NRDC  database.  The  sample  rate  for 
these  data  is  250sps,  and  many  events  are  believed  to  be  commercial  explosions.  In  one 
example,  a  particular  linear  filter  was  found  which  transforms  the  data  segment  (Figure  la) 
into  a  series  of  irregularly  spaced  impulses  (Figure  lb).  These  impulses  have  interphase 
arrival  delays  of  the  correct  order  of  magnitude  for  regional  phases  at  this  distance.  Note 
that  the  time  separations  between  the  doublets  are  uniform  (Figure  2). 

The  data  analyzed  from  the  Soviet/NRDC  database  were  recorded  with  high-gain  sur¬ 
face  instruments.  Figure  la  shows  data  recorded  at  station  Karkaralinsk  (KK)  for  an  event 
at  (49.9N,  73.  IE),  at  a  distance  of  1.6  degrees.  The  data  window  shown  is  approximately 
6.0  seconds  of  data,  starting  just  before  the  Pg  arrival.  A  2.8-second  (700  point)  MED 
filter  was  generated  for  this  data  window.  The  filtered  data  is  shown  in  Figure  lb.  The 
MED  algorithm  does  not  constrain  the  phase  of  the  filter,  and  the  output  of  the  filter 
operation  is  therefore  time  shifted  by  an  arbitrary  amount.  The  filtered  data  contains  four 
dominant  spikes,  with  time  separations  of  .34  sec,  .16  sec  and  .52  sec.  Three  of  the  spikes 
appear  as  doublets,  uniformly  separated  by  47  msec,  suggesting  source  multiplicity. 

We  obtained  data  at  the  end  of  the  second  quarter  from  a  recent  known  quarry  blast 
in  Oklahoma,  along  with  a  description  of  the  shot  geometry  and  delay  times.  The  shot 
description  is  shown  in  Figure  3.  The  typical  shot  delay  time  is  approximately  8msec. 
Since  the  data  sample  rate  is  60sps  (16.67msec  sample  interval),  we  can  not  deconvolve 
the  individual  shots  in  this  case.  Figure  4  shows  a  spectrogram  of  the  first  1.5  seconds  of 
a  signal  from  the  quarry  blast,  recorded  with  a  borehole  instrument  approximately  32km 
from  the  quarry.  There  are  two  prominent  spectral  lines  at  approximately  42msec 
(23.8Hz)  and  59msec  (16.9hz),  which  would  appear  to  be  directly  related  to  the  shot  delay 
times. 

Figure  6  shows  the  effect  of  a  50  sample  MED  filter  on  the  first  1.2  seconds  of  the 
signal,  as  shown  in  Figure  5.  The  deconvolved  signal  contains  4  spikes  with  time 


-  1  - 


Figure  2.  The  deconvolved  signal  from  Figure  lb  enlarged  to  show  the  doublet  nature  of 
the  impulses. 
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Figure  3.  Quarry  explosion  timing  and  layout. 
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Figure  4.  Spectrogram  of  explosion  example  showing 
peaks  at  16.6  and  23.8Hz  discussed  in  text. 


input  data:  LNO/szl  tiller  output 


Figure  5.  Explosion  signal  used  as  input  to  Med  filter  Figure  6.  Output  from  50-sample  MED  filteron  first  1 .2 

seconds  of  signal  shown  in  Figure  3. 
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separations  of  52  msec,  1 19  msec,  and  580  msec.  This  result  may  be  related  to  the  source 
multiplicity. 

3.  FUTURE  PLANS 

Further  analysis  of  the  data  sets  with  which  we  are  currently  working  should  deter¬ 
mine  whether  any  of  the  impulse  signals  we  obtain  with  the  MED  filter  are  related  to 
regional  phase  arrivals.  We  intend  to  acquire  data  from  known  quarry  explosions,  sam¬ 
pled  at  a  higher  rate  than  the  data  with  which  we  are  currently  working,  to  test  further  the 
usefulness  of  this  method  in  detecting  source  multiplicity. 
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1 

AWRE 

MARSHALL 

i 

NOR 

UARTN 

i 

FRANCE.  RADICMANA 

UASSNON,  kCCHBEH 

2 

SMWHYSCSOST 

MCCARTOHGRAY 

1 

8-CUBED,  LA  JOLLA 

MCLMAGKN 

1 

UCBD 

MNSim  ORCLITT.  GIVEN 

2 

ST  LOUIS  U 

MITCHBi,  HERRMAAN 

1 

SCUBED.RESTON 

MLRPW 

2 

RADIX 

PULLI 

1 

NORWAY.  NTNF 

nnggal 

2 

TEL  EDYNE.  ALEXANDRIA.  VA 

RIVERS 

2 

TASC 

SAILOR 

1 

MAY-2 1-1S91  08=25  FROM  DARPA 


TO 


91505846960?  P.03 


CO fit  UAXJNG  UST-MM 


CROAHttATON 

HUC 

NO.  COPIES 

use 

SAUMIS.  AKI 

1 

IRIS 

SMP90N 

2 

U.  WASKNQTON 

SMITH 

1 

UWSCCNSM 

THJRFSSUCYST 

1 

MIT 

TOK302A3AWTY 

1 

U.AZ 

WALLACE 

1 

MRC.NEWNGTCHYA 

WCRTMAN 

" 

1 

USQOVtRJrtCXrAOtX** 

ACOA 

UBBERMAN 

1 

ATOSRNP 

jBwrcmcp 

1 

AFTAC.  CSS.  R06SLYN.  VA 

blandtcrd 

1 

AFTAC/CA 

S7TJPO 

1 

AFTAC/TT 

PLOTTE 

1 

ClA/ACtS 

KATIE  FOLEY 

1 

ClAOSWR 

TUfUaiJL 

1 

OAflPA 

ALEWWE  RYALL  KERR 

7 

OARPAFWO 

LIBRARIAN 

1 

01A 

GLCV53 

1 

0NA/8PSS 

9CRE 

1 

tee 

KDCNTZ 

1 

one 

NPOCTR 

Z 

GULWH 

CIPAR 

1 

GU.WH 

LEWMowicz 

1 

<317300 

30 

1 

LLNL 

HYSNCN 

2 

ago 

CERE 

1 

SANOtA 

CHAR 

1 

IBS 

LBTH 

1 

U333 

MASSE 

1 

WUKTBSO 

FENOE 

1 

TOTAL  NUMBER  OF  REPORTS 

88 

V 
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